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-^A  conical  shell  finite  element  is  used  to  model  the  shell 
and  the  exterior  Helmholtz  integral  equation  the  fluid,  the 
displacements  and  pressure  being  represented  by  Fourier  series 
expansions  in  the  circumferential  coordinate.  The  dynamic  stiff¬ 
ness  matrix  of  the  system  is  obtained  by  combining  the  separate 
dynamic  stiffness  matrices  of  the  fluid  and  structure,  using 
'normal'  coordinates  alone  to  save  computer  time.  The  excitation 
is  either  a  time-harmonic  point  force  or  a  plane  wave  incident  at 
an  oblique  angle.  In  order  to  validate  the  problem  formulation 
and  computer  programs,  numerical  results  of  far-field  sound 
radiation  from  a  spherical  elastic  shell  are  compared jwith  those 
from  a  closed-form  shell  theory.  C  •C'-?  "%  ' •  - 
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INTRODUCTION 


A  previous  Technical  Memorandum  C1J  describes  the  method  of 
dynamic  stiffness  coupling,  using  normal  coordinates  alone  to 
couple  the  matrix  equations  of  motion  of  the  structure  and  fluid, 
as  a  potentially  useful  tool  for  solving  radiation  and  scattering 
problems  of  submerged  axisymmetric  shell  structures  which  are 
subjected  to  axisymmetric  excitation.  The  conical  shell  element 
was  used  model  the  shell  and  the  interior  Helmholtz  integral 
equation  the  exterior  fluid.  It  recommended  that  the  work  should 
be  extended  to  include  the  case  of  general  excitation  in  which  the 
field  quantities  of  interest  are  expanded  as  Fourier  series  in  the 
the  circumferential  coordinate. 

It  is  the  aim  of  the  work  contained  herein  to  present  the 
mathematical  formulation  of  the  problem  of  the  acoustics  of 
axisymmetric  shell  structures,  which  are  subjected  to  asymmetric 
excitations,  and  to  demonstrate  numerically,  by  comparisons  with 
an  exact  theory,  that  a  particular  Fortran  program  implementation 
is  a  useful  tool  for  research  studies. 

For  completeness,  this  Technical  Memorandum  includes  much 
of  the  work  contained  in  the  previous  study  C13,  which  public¬ 
ation  is  now  out  of  date.  The  derivation  of  the  matrix  equations 
of  motion  is  given  without  elaboration  and  with  minimal  explan¬ 
ation,  and  a  good  working  knowledge  of  the  finite  element  method 
and  the  Helmholtz  integral  equation  is  required  for  a  proper 
understanding  of  the  Memorandum's  contents. 


2.  FINITE  ELEMENT  METHOD 


A  finite  element  Fortran  program  C23  i3  available  for  the 
dynamic  analysis  of  axisymmetric  thin  shell  structures  whose 
displacements  are  represented  by  Fourier  series  in  the  circumfer¬ 
ential  coordinate,  viz. 

u(r,z,$)  =  r°°  „  u  (r,z)cos(n$) 
n=u  n 


v(r,z,4>) 


r“  v  (r,z)sin(n<p) 
n=u  n 


V(r,z,<ti) 


.  w  <  r,z)cos(n$) 
n=u  n 


(2.1) 


5<r,z,<>)  =  z"=0  5n(r,z)cos(n$) 

in  which  u,  v,  w  are  axial,  tangential  and  radial  displacements, 
and  £  is  the  meridlanal  rotation.  Time  variation  exp(-iu)t)  is 
to  be  omitted  throughout. 


The  conical  shell  element  is  the  basic  structural 
component,  its  geometry  being  shown  in  Figure  1c.  The  mass  and 
stiffness  matrix  in  the  n'th  circumferential  harmonic  have  been 
given  by  Percy  et  al  C33.  Assembly  of  the  elements  gives  the 
system  matrix  equation 
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(2.2) 


CS  -w2M  3  CX  3  =  CF_3 
n  n  n  n 

4m*4m  4m*l  4m*l 

for  each  harmonic,  where  m  Is  the  number  of  nodes  on  the 
structure;  CSr3  and  CMn3  are  the  s7stein  stiffness  and  mass 

matrices.  Numerical  inversion  of  this  equation  is  rapid  if  the 
bandwidth  is  kept  small.  Thus,  given  the  nodal  force  excitation 
vector  CFn3,  the  response  vector  CX^  can  usually  be  found 

without  excessive  computer  run-times. 

The  Fortran  program  has  been  modified  to  give  the  'normal' 
receptance  matrix,  which  is  defined  as  the  normal  displacement  at 
surface  point  'i'  due  to  unit  'normal'  force  at  surface  point 
'j'.  If  CW3  and  CF3  are  now  defined  as  matrices  of  'normal' 
displacements  and  'normal'  forces,  in  the  n'th  circumferential 
harmonic,  at  the  m,  say,  surface  nodes,  then 

CW3  =  CS  (U))3'1CF3  (2.3) 

e 

ra*l  m*m  m*l 

where  CSe(w)3  1  is  the  frequency  dependent  receptance  matrix. 

The  inverse  of  this  equation  is  the  dynamic  stiffness  relation 

CSe(w)3  CW3  =  CF3  (2.4) 

m*m  m*l  m*l 

which  is  the  starting  point  of  the  fluid-structure  interation 
matrix  analysis. 


3-  HELMHOLTZ  INTEGRAL  EQUATION 
(a)  General 


The  Helmholtz  integral  equation  for  a  smooth  surface  Sg 
undergoing  time-harmonic  oscillations  in  an  acoustic  fluid  is 

q(R)p(R)  =  (l/4ir)  /  Cp(  RQ  )  3G(  B  ,E0 )  /  3NQ  -  pw2W(  8Q  )G<B,B0  HdSg 

(3.1) 

where  p  is  acoustic  pressure  and  W  is  normal  surface  displace¬ 
ment.  q(R>  is  1.0,  0.5  or  0  depending  on  whether  R  is  exterior 
to,  on,  or  Inside  the  closed  surface  Sq.  3/3Nq  is  different¬ 
iation  in  the  direction  of  the  outward  normal. 

G(R,R0>  -  exp(ik|R-R0l )/|R-R0l  (3.2) 

is  the  free-space  Green's  function;  p  is  the  density  of  the 
exterior  fluid  whose  sound  velocity  is  c  and  whose  acoustic 
wavenumber  is  k»  w/c.  The  geometry  is  shown  in  Figure  la. 
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Problems  associated  with  a  proper  implementation,  which 
is  unique  at  all  frequencies,  have  been  analysed  by  Burton  for 
the  separate  cases  of  general  C43  and  axisymmetric  surfaces  C53. 
Computer  implementation  is  not  straightforward  due  to  the  presence 
of  strongly  singular  integrals,  which  are  difficult  to  evaluate 
numerically.  It  is  the  author's  impression  that  implementations 
discussed  in  the  literature  either  ignore  or  do  not  treat  properly 
these  singularities.  Because  the  integrands  are  integrable  the 
singularities  may  be  ignored,  but  not  safely  because  the  resulting 
computer  programs  may  be  unreliable  as  well  as  inefficient.  A 
recent  publication  by  Reut  C63  presents  formulae  which  remove  the 
problem  of  singular  integrals  for  the  case  of  flat  rectangular 
surface  elements  over  which  the  pressure  and  displacement  can  be 
regarded  as  constant. 

(b)  Axi symmetric  Surface 

For  the  specific  case  of  the  axisymmetric  surface  with  the 
pressure  p  and  normal  displacement  W  represented  by  Fourier 
series  expansions  in  the  circumferential  coordinate,  viz. 

p(r,z,0)  =  s“_0  pn(r,z)cos(n$) 

(3.3) 

W(r,z,<f>)  =  Wn<r ,z)cos(n$) 

a  procedure  which  explicitly  takes  into  account  the  nature  of  the 
singularity  in  the  integrands  has  been  recommended  by  Miller  C73. 


The  Green's  function  is  expanded  as 
G(R,B0>  =  2"=0  Gn(r,r0,z-z0,)cosn(<>-<|>0>  (3.4) 

in  which 

G  (r,r„,z-z„)  =  (e  /*)  /*  Cexp(  ikD)  /D3cos(n<j>>d<|>  (3.5) 

n  u  o  n  u 

with  en=l  for  n=0,  en=2  for  n*l,2,3,..;  and 

D  =  D(r,r0,z-zQ)  =  £( z-z0)2+r2+rg-2rr0cos($) 31/2  (3.6) 

The  normal  derivative  of  the  Green's  function,  viz. 

9G/dN0  -  E"a0  C3Gn/3No3cosn(0-0o)  (3.7) 

is  obtained  from  the  formula 

3G  /3N.  *  cos( B) 3G  /3r_  -  sin< P ) 3G  / 3zn  (3.B) 

no  no  n  u 


in  which  the  'normal'  angle  p  is  defined  in  Figure  la.  Thus, 
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in  -snsa 


3Gn/3N0  =  (en/ir>  /q  <i»c-l/D)Cexp(iW3)/D2}  X 

C (r0~r .cos$)cosp-(zo-z>sinp}cos(n$)d0  (3.9) 

Noting  that  the  area  element  on  the  surface,  after  first 
performing  the  integration,  is 

dSQ  =  (2ir/en)  (rQ/cosP)dz0  (3.10) 

enables  the  Helmholtz  integral  equation  to  be  written  in  terms 
of  the  harmonic  amplitudes  as 

q(r,z)pn(r,z)  =  <l/2en)/  Cpn< rQ , zQ ) 3Gn< r ,rQ ,z-zQ ) / 3NQ 

-pw2W(r0,z0)Gn(r,r0,z-z0) J(r0/cosp)dz0  (3.11) 

The  integral  runs  over  the  axial  length  of  the  surface. 

( c )  Discretization 

Let  the  field  point  (r,z)  be  on  the  surface,  which  is 
divided  into  m-1  bands,  as  in  Figure  lb,  with  ends  at  z^,  i=l,m. 

Thus  q(r,z)  is  0.5  in  equation  (3.11),  giving 

m-1 

-Pn(r,z)  +  <l/en)£  /  Cpn(r0,z0)3Gn(r,r0,z-z0)/3N0J(r0/cosP)dz0 
i=l 

m-1  (3.12) 

=  (pw2/en)Z  I  W(r0,z0)Gn(r,r0,z-z0>(r0/cosp)dz0 
i  =  l 

Along  the  generator  of  any  element  let  the  acoustic  pressure 
and  normal  displacement  vary  linearly  with  the  axial  coordinate, 
zQ,  viz. 


Pn(r0'V  =<zi+l  2 


Vr0'V  =(zi+l"zi)  Czi+rv 


i}  lczi+l  z0'-zi+z0^  Tpi  1 

Lpi+lJ 

iJ  1^zi+l-z0,-zi+z03  j*Wi  j 


(3.13) 


where  p.^  and  pi+1  are  the  values  of  the  surface  pressure 
harmonic,  p^,  at  the  ends  of  the  element;  and  and  are 

the  values  of  the  normal  displacement  harmonic,  Wn,  at  the 
ends  of  the  element. 


Substitute  equation  (3.13)  into  equation  (3.12)  to  give  the 
matrix  equation 


m-1 

-pn<r,z)+  Z  CCa(r,z) ,b(r,z)]|pi  1-  Cc(r,z) ,d(r,z)3rw. 

1=1  kJ  LH 


i 

i+1 


(3.14) 

In  which 

a(r,z)  =  h  /  (zi+1-z0)£3Gn(r,r0,z-z0)/3N0l(r0/cosP)dz0 


b(r,z)  =  h  /  <-z.+zn)£3G  (r.r.fZ-z.WSNrtKr./cc^Jdz 


c(r,z)  =  pw  h  /  (zi+1-z0)Gn(r,r0,z-zQ) (rQ/cosp)dzc 


d(r,z)  =  pu  h  /  (-zi+z0)Gn(r,rQ,z-z0> (rQ/cos3)dz0 


(3.15) 


,-1 


with  h=(zi+1-zi>  /en,  the  integration  limits  being  z ^  and  zi+1- 


If  equation  (3.14)  is  written  down  at  the  m  values  of  (r,z) 
which  correspond  to  the  surface  nodes,  then  a  matrix  equation 
results  which  relates  the  surface  pressure  and  displacement  at  the 
m  surface  nodes,  viz. 


CB(u>)  3  Cp3  =  CC(w)3  CW3  (3.16) 

m*m  m*l  m*m  m*l 

The  inverse 

Cp3  =  CB(w)3~1CC(w>3  CW3  =  CD(w>3  CW3  (3.17) 

m*l  m*m  m*m  m*l  m*m  m*l 


of  this  equation  gives  the  surface  pressure  matrix  in  terms  of  the 
normal  displacement  matrix. 

The  integrals  in  equation  (3.12)  must  be  evaluated  with  some 
care  as  they  are  singular  (but  integrable).  Details  of  their 
numerical  evaluation,  based  on  Miller's  work  C73,  are  to  be  given 
elsewhere. 


4.  COUPLED  FE  AND  HIE 

The  matrix  equation  of  motion  is 

CS  <w)3  CW3  *  CF-3  +  CF„3  (4.1) 

e  E  p 

m*m  m*l  m*l  ra*l 

where  CF-.3  are  the  mechanical  excitations,  and  CF  3  are  the  forces 
E  p 

due  to  fluid  pressure  which  are  obtained  from  the  nodal  point 
fluid  pressures  by  the  relation 
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(4.2) 


CFp3  =  -CA3  Cp3 
m*l  m^m  m*l 

where  CAD  is  a  frequency  independent  area  matrix  that  converts 
nodal  pressures  to  nodal  forces. 

The  equivalent  nodal  point  forces,  due  to  fluid  pressure, 
on  the  ith  element  are  obtained  by  a  standard  finite  element 
procedure.  For  linear  variation  of  both  displacement  and 
pressure  they  are 

=  -<2ir/en>  /  c3  cxc2  p±  1  ( rQ/cosp>dz0 

_clc2  c2  _  fi+1  <4*3> 

where 


C1  *  <zi+l-zi>  (2i+i-2o> 


c2  =  ^i+l'V  (-zi+20> 


(4.4) 


and  the  range  of  integration  is  2 ^  to  Assembly  of  the 

element  force  matrix,  equation  (4.3),  over  all  the  elements  gives 
the  frequency  independent  triple-diagonal  area  matrix,  CAD. 

Combining  equations  (4.1)  and  (4.2)  gives  the  matrix 
equation  of  motion  as 

CSe<w)D  CW3  *  CFe3  -  CAD  CpD  (4.5) 

m*m  m*l  m*l  m*m  m*l 


which  is  now  made  specific  to  either  mechanical  excitation  or 
acoustic  excitation. 


5.  MECHANICAL  EXCITATION 
(a)  Response 

Expand  the  mechanical  excitation  stress  in  the  normal 
direction  as  a  Fourier  series,  viz. 

F(r0,z0,#0)  =  2"_0  Fn(r0,z0)cos(n$0)  (5.1) 

where 

Fn(r0'z0)  *  (en/2w)  SV  F<ro,zo,Oo)cos(n0o)d(|>o  (5.2) 

For  the  specific  case  of  a  point  force,  FQ,  located  at  one  of 
the  nodes,  j  say,  at  4*0 

F<r0,z0,*0>  •  F0$($>6{z0-z^)(cosp/r0>  (5.3) 
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the  harmonic  amplitudes  are  obtained  from  equation  (5.2)  as 


Fn(r0,z0>  =  <enFocos0/2irro>S(  zQ-z ^  )  (5.4) 

which  can  be  shown  to  correspond,  using  equation  (4.3)  and  the 
standard  sign  convention  which  requires  positive  values  of 
stress  and  pressure  to  act  in  opposite  directions  at  a  surface, 
to  a  force  vector,  CFg3,  which  has  all  zero  elements  except 

for  the  j'th  which  is  simply  Fq  for  all  harmonics. 

The  fluid  pressure  caused  by  the  shell's  vibrations  must 
satisfy  the  the  Helmholtz  integral  equation;  and  thus  the  matrix 
relation  equation  (3.17),  which  when  substituted  into  equation 
(4.5)  gives  the  coupled  equations  of  motion  as 

|fSe(w)3+CADCD(w)DjcW3  =  CF£3  (5.5) 

m*m  m*l  m*l 

and  it  should  be  noted  that  the  coefficients  of  the  m*m  matrix 
are  neither  symmetric  nor  banded,  because  the  matrix  CD(w)3  has 
been  derived  from  an  integral  equation  rather  than  an  energy 
variational  principle. 

Inversion  of  the  matrix  equation  (5.5)  gives  the  'normal' 
displacement  response  CW3  in  the  n'th  harmonic.  The  surface 
pressure  Cp3  in  the  n'th  harmonic  is  found  from  equation  (3.17). 

<b)  Far-Fleld  Sound  Radiation 


Approximating  the  Green's  function  and  its  derivative, 
in  equations  (3.5)  and  (3.9)  by  their  values  for  a  large  'D' 
gives  the  stationary  phase  approximation  to  the  far-field 
pressure  in  the  n'th  harmonic  (after  much  algebra)  as 

pnf(R,0,<t>>  =  Y^a)  (-Dnexp(ikR)/2R  (5.6) 

where 

Y(a) 

with 

Y(zQ>  =  (r0/cosp)  X  (5.8) 

[pn<r0'20)Cv-cosP'Jn(^r0>+i(x'sinPJn<'Yr0>:>'pw2wntr0'20)Jn(Yr0)] 

-Y=lc.sin9  and  <x=)c.cos0  being  the  stationary  phase  wavenumbers. 

The  integral  in  equation  (5.7)  is  lasily  evaluated  by,  for 
example,  Gaussian  quadrature,  the  surface  pressure  pn  and  the 

surface  normal  displacement  Wn  being  found  by  interpolation  in 

a  table  of  Cp3  and  CW3  values  which  are  found  as  described  above. 


/  Y(z0)exp(-i»z0)dz0 
2 1 


(5.7) 


t. 


6.  ACOUSTIC  EXCITATION 
(a)  General  Excitation 

It  is  convenient  to  express  the  total  fluid  pressure  in 
the  n'th  harmonic  as 

Pn<r,z)  =  Plnc<r,z)+pr<r,z)+pe<r,z)  (6.1) 

where  Plnc  is  the  incident  wave,  p  is  the  acoustic  pressure 
scattered  as  though  the  shell  were  rigid  (and  immovable)  and  pg 

is  the  pressure  scattered  by  the  shell's  elastic  motion.  The 
total  acoustic  particle  displacement,  Wn(r,z>,  in  a  direction 

normal  to  the  surface  is  found  from  the  acoustic  momentum 
relation 

2 

3pn(r,z)/3N0  =  cosp. 3pn< r , z ) /3r-sinS. 3pn ( r ,z ) / 3z  =  pu>Wn<r,z) 

(6.2) 

The  nodal  pressures  and  particle  displacements  are  given  by 
the  matrix  relations 


CP3  =  Cp.nc3+Cprl+Cpe3 
m*l  m*l  m*l  m*l 

CW3  =  CW.  3+CW  3+CW  3  =  CW  3 
inc  re  e 

m*l  m*l  m*l  m*l  m*l 


(6.3) 


Note  that  the  rigid  body  pressure  is  such  that  CWr3  =  ~CWinc3, 

which  means  that  it  can  be  thought  of  as  that  pressure  which 
arises  from  a  virtual  displacement  -CWinc3  normal  to  the 

surface.  The  particle  displacement  CWin(,3  is  found  from  the 
discretized  version  of  equation  (6.2)  in  which  pn<r,z)  is  set 
to  Pinc<r,z)  and  W  ( r ,  z )  to  Win(.(r,z).  Thus  CPinc3r  CWinc'J 
and  hence  CWr3  are  all  known. 


The  nodal  surface  pressure  and  displacement  of  the  rigid 
body  motion  (suffix  r)  and  the  elastic  motion  (suffix  e)  must 
be  solutions  of  the  Helmholtz  integral  equation.  Thus  from 
equation  (3.17) 

Cpr3  =  CD(w>3  CWr3  (6.4) 

*1*1  m*m  m*l 

is  the  (now)  known  matrix  of  of  'rigid'  nodal  pressures,  and 

Cpe3  =  CD(w)3  CWe3  (6.5) 

m*l  m*m  m*l 
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is  the  matrix  relation  which  connects  the  unknown  'elastic' 
nodal  pressures  and  displacements. 

From  equation  (4.5),  setting  CFED=0,  together  with  the 
first  of  equations  (6.3) 

CSe(w)D  CWD  =  -CAD  cPlnc+Pr+Pe3  (6.6) 

Substitute  equation  (6.5)  into  this  equation  to  give  the 
fundamental  fluid-structure  matrix  equation  of  motion  for 
sound  scattering  as 

[cSe<w]|+CADCD(u))D}cW3  =  -CAD  Cpinc+Pr3  (6.7) 

m*m  m*l  m*m  m*l 


and  it  is  evident  that  for  the  case  of  acoustic  excitation 
the  driving  force  is  the  equivalent  force  due  to  the  'blocked' 
acoustic  pressure. 

Inversion  of  equation  (6.7)  gives  the  total  nodal  point 
displacement  CWD  =  CWeD,  and  hence  the  'elastic'  nodal  pressure 

matrix,  CpeD,  from  equation  (6.5).  The  scattered  pressure 

and  its  normal  displacement  at  the  surface  nodes  are  simply 
defined  as 


Cpg]  =  CprD+CpeD 

CW  D  =  CW  D+CW  D  =  -CW.  D+CKID 
s  re  me  e 


(6.8) 


The  pressure  scattered  to  the  acoustic  far-field  is  found  from 
equations  (5. 6-5. 8)  in  which  Pn(ro'zo'  and  Wn<r0'z0)  are  found 

by  interpolation  in  the  table  of  nodal  point  values  of  equation 
(6.8).  A  frequently  used  reference  is  the  pressure  scattered  to 
the  far-field  by  the  shell  assumed  to  be  rigid:  this  is  obtained 
from  equations  15.6-5.8),  using  Cp  D  and  CW  3. 


(b)  Plane  Wave  Excitation 


Let  the  shell  be  insonified  by  the  plane  wave 
pinc<x'y,z)  =  Pj^expfikfx.  sinGj-z.cose^  )  }  (6.9) 

which  is  incident  at  the  angle  $^=180°,  see  Figure  la.  Its 
expansion  in  cylindrical  coordinates  is  C8D 


Pinc<r,z,$)  *  Pjexpl -ikz . cos©^ )l“_0  eninJn( kr . sinQj ) cos (n$ ) 

(6.10) 

From  now  on  the  summation  sign  and  the  cos(n$)  factor 
will  be  dropped  because  the  analysis  herein  refers  to  the 
harmonic  amplitudes  only.  Thus  the  incident  pressure  is 
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(6.11) 


pinc(r,z>  =  piexP* -ikz. cos©^ )eninJn(kr . sin6^ ) 

whose  particle  displacement  is  obtained  from  equation  (6.2), 
with  Pinc  substituted  for  pn,  as 

Winc(r'z)  =  < l/pw2)exp( -ikz. cos6i )enin  X  (6.12) 

I 

£k.  sin0^Jn<  kr .  sin6^ )  cosp+ik.  cos6^ ( kr .  sine^^ )  sin61 

All  the  ingredients  are  now  present  to  enable  computation 
of  the  scattered  pressure  field.  The  procedure  to  be  followed 
is  the  same  as  for  the  general  excitation  case,  which  is 
described  above. 


7 .  NUMERICAL  TESTS 
(a)  General 


Fortran  programs  have  been  written  for  calculating  the 
far-field  pressure  from  an  axisymmetric  shell  structure  (with 
or  without  branches  and  ribs)  which  is  excited  by  either  a  time- 
harmonic  mechanical  point  force,  or  by  a  time-harmonic  plane  wave 
which  is  incident  at  an  arbitrary  elevation  angle.  The  field 
quantities,  of  course,  are  represented  by  Fourier  series  expan¬ 
sions  in  the  circumferential  coordinate. 

Miller  C7]  has  shown  that  that  the  integral  represent¬ 
ations  of  the  Green's  function  and  its  derivative,  equations 
(3.5)  and  (3.9)  introduce  logarithmic  singularities  in  the  kernel 
of  the  integral  equation  (3.11),  at  the  point  (r0,zQ)=(r,z) ,  and 

goes  on  to  develop  quadrature  formulae  which  explicitly  treat 
these  singularities.  However,  because  the  singularities  are 
logarithmic,  they  may  be  integrated,  although  inefficiently, 
by  Gaussian  quadrature.  The  penalty  of  course  is  that  the 
computation  time  needed  to  achieve  a  given  accuracy  increases 
several  tines.  The  current  versions  of  the  programs  use 
Gaussian  quadrature,  thus  effectively  ignoring  the  singularity 
except  for  vastly  increasing  the  quadrature  points  on  singular 
elements.  A  future  publication  will  show  the  value  of  treating 
these  singularities  properly  and  it  will  also  describe  an 
attempt  to  improve  the  integration  over  end  elements,  which 
turned  out  to  be  troublesome. 

(b)  Spherical  Shell  Constants 

As  a  first  test  of  the  correctness  of  the  Fortran  programs 
numerical  results  of  sound  radiation  from  and  sound  scattering  by 
a  thin  empty  spherical  steel  shell,  submerged  in  water,  are 
compared  with  those  obtained  from  closed-form  expressions  derived 
from  a  shell  theory  C9,103.  The  material  and  geometric  constants 
in  SI  units  are  as  follows: 

Shell:  Young's  modulus  19.5x10^® 

Poisson's  ratio  0.29 
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Density  7700.0 
Thickness  0.01 
Radius  1.0 

Hysteretic  loss-factor  0.01 

Water:  Density  1000.0 

Sound  velocity  1500.0 

< c )  Spherical  Shell  Test  of  n=0 

Figure  2  shows  radiated  sound  levels  when  the  shell  is 
excitated  by  a  'normal'  point  force  located  on  the  z-axis. 

Thus  both  the  geometry  and  the  excitation  are  axisymmetric , 
which  means  that  only  the  n=0  harmonic  contributes  to  the  far- 
field  radiation,  calculated  at  6=0.  In  Figure  2a  the  shell 
was  divided  into  20  conical  shell  elements,  with  a  constant 
angular  spacing  of  9°.  Agreement  with  the  'exact'  solution  is 
good  up  to  the  frequency  of  the  first  resonance;  thereafter  the 
finite  element  levels  are  seriously  in  error,  although  the 
frequencies  at  which  resonances  occur  are  reasonably  accurate. 

In  Figure  2b  the  shell  was  divided  into  40  conical  shell 
elements  with  an  angular  spacing  of  4.5°.  Agreement  with  the 
exact  solution  has  improved  dramatically. 

Figure  3  shows  plots  of  monostatic  target  strength,  the 
excitation  being  an  axial  plane  wave.  Again  both  the  geometry 
and  the  excitation  are  axisymmetric,  which  means  that  only  the 
n=0  harmonic  contributes  to  the  target  strength.  In  Figure  3a 
the  shell  was  divided  into  20  conical  shell  elements,  with  a 
constant  angular  spacing  of  9®.  Agreement  with  the  'exact' 
solution  is  very  good  up  to  the  frequency  of  the  second  reson¬ 
ance;  thereafter  the  finite  element  levels  are  in  error,  but  not 
as  seriously  as  those  in  Figure  2a.  The  shell  was  divided  into 
40  conical  shell  elements  with  an  angular  spacing  of  4.5®  in 
Figure  3b.  Agreement  with  the  exact  solution  has  improved,  as 
expected,  and  is  satisfactory  over  all  the  frequency  range. 

Also  shown  in  Figure  3  is  the  sound  scattered  by  a  rigid 
sphere.  There  is  no  distinguishable  difference  between  the 
'exact'  curve  and  the  levels  derived  (in  this  case)  from  the 
Helmholtz  integral  equation  alone. 

The  numerical  results  discussed  have  demonstrated  that  the 
problem  formulation  and  computer  programs  must  be  substantially 
correct,  at  least  for  the  n=0  harmonic.  They  also  suggest  that 
the  weak  link  in  the  numerical  formulation  is  more  likely  to  be 
in  the  finite  element  part  rather  than  in  the  Helmholtz  integral 
equation  implementation.  No  substantiated  explanation  is  offered 
as  to  why  the  numerical  results  of  sound  scattering  are  more 
accurate  than  those  of  acoustic  radiation,  for  a  given  idealis¬ 
ation.  Perhaps  it  is  related  to  the  accuracy  of  the  off-diagonal 
elements  of  the  shell's  receptance  matrix  being  more  accurate 
than  the  diagonal  elements. 

<d)  Spherical  Shell  Test  of  n=0-3 

The  implementation  of  the  Helmholtz  integral  equation  should 
fail  at  the  forbidden  frequencies,  which  are  the  eigenvalues  of 
the  associated  interior  problem  with  Dirichlet  boundary  condition 
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on  the  surface.  Practically,  this  means  that  the  numerical 
inversion  of  the  matrix  CB3  defined  by  equation  (3.16)  will 
become  more  and  more  ill-conditioned  as  a  forbidden  frequency  is 
approached.  For  the  sphere  the  forbidden  frequencies  are  simply 
the  roots  of  the  equation,  jn<ka)=0,  which  are  f=750  for  n=0, 

f=1072  for  n=l,  f=1376  for  n=2  ...  Numerical  tests  did  indeed 
show  severe  111-conditioning  as  these  critical  frequencies  were 
approached . 

As  a  test  of  the  full  formulation  for  general  'n'  it  is 
useful  to  consider  the  sound  scattered  by  the  spherical  shell 
when  it  is  excited  by  a  plane  wave  incident  at  0^90°.  In  this 

case  the  monostatic  target  strength  is  a  combination  of  all 
harmonics.  Because  of  symmetry,  it  must  be  the  same  (at  least 
theoretically)  as  the  monostatic  target  strength  of  the  shell 
when  it  is  excited  at  0=0°,  which  involves  the  n=0  harmonic 
alone.  The  following  table,  in  which  R  denotes  scattering  from 
a  rigid  body  and  E  denotes  scattering  from  the  spherical 
elastic  shell,  is  a  comparison  of  monostatic  target  strength 
levels  in  decibels  at  selected  frequencies: 


Frequency  Hz 

Exact 

R  E 

CD 

II 

»o 

o 

(n=0) 

E 

0  =  90° 
R 

(n=0-3 ) 
E 

100 

-17.5 

-19.9 

-17.5 

-19.9 

-17.5 

-19.9 

150 

-11.5 

-13.0 

-11.5 

-13.0 

-11.5 

-13.0 

200 

-8.0 

-7.6 

-8.0 

-7.7 

-8.0 

-7.6 

250 

-6.3 

1.0 

-6.4 

0.4 

-6.3 

0.9 

300 

-6.2 

-4.8 

-6.3 

-3.8 

-6.2 

-4.2 

For  6=90°  numerical  tests  showed  that  it  was  sufficient 
to  sum  from  n=0  to  3  only.  The  main  contribution  comes  from 
the  n=l  harmonic  and  the  n=0  harmonic  is  significant  also;  the 
n=2  harmonic  contributes  significantly  at  200,  250  and  300Hz, 
and  the  n=3  harmonic  is  significant  only  at  300Hz.  The  'rigid' 
body  levels  (R)  agree  to  within  O.ldB  at  all  frequencies,  which 
demonstrates  that  the  solution  of  the  Helmholtz  integral  equation 
is  also  credible  at  values  of  'n'  other  than  n=0.  The  spherical 
shell  levels  (E)  agree  to  within  O.ldB  at  100,  150  and  200Hz,  and 
to  within  ldB  at  frequencies  250  and  300Hz  which  are  either  side 
of  a  resonance.  Overall  these  results  help  considerably  to  rein¬ 
force  the  opinion  that  the  programs  could  be  a  useful  tool  for 
solving  certain  low  frequency  radiation  and  scattering  problems. 


8.  CONCLUDING  REMARKS 

A  receptance  method  for  calculating  the  exterior  acoustics 
of  axisymmetric  shell  structures  has  been  presented.  The  method 
is  relatively  economical  in  computer  time  because  'normal'  coord¬ 
inates  alone  are  used  for  the  fluid-structure  coupling.  The 
method  has  been  validated  by  comparing  numerical  results  obtained 
from  the  formulation  with  'exact'  results  obtained  from  a 
closed-form  expression. 

Much  theoretical  work  remains  to  be  done  before  a  practical 
and  robust  suite  of  programs  is  available,  but  it  is  thought  that 
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the  programs  developed  so  far  are  a  useful  tool  for  obtaining 
numerical  results  of  sound  radiation  and  scattering  for  research 
purposes.  Future  work  in  this  area  should  include  (a)  a  proper 
treatment  of  the  singular  integrals  in  the  Helmholtz  integral 
equation,  according  to  Miller's  analysis  C73;  and  an  implementat¬ 
ion  of  a  quadratic  isoparametric  element  to  replace  the  linear 
element  used  herein;  (b>  a  theoretical  analysis  of  the  additional 
singular  integral  which  would  occur  if  it  were  thought  desirable 
to  Implement  the  full  Burton  and  Miller  method,  described  in 
Burton  C4,5J;  and  (c)  the  use  of  a  quadratic  isoparametric  conical 
shell  element,  as  the  numerical  results  have  suggested  that  the 
finite  element  part  is  the  weakest  link  in  the  formulation. 
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